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Abstract 

A search for pair-production of neutralinos at a LEP centre-of-mass energy 
of 189 GeV gave no evidence for a signal. This limits the neutralino produc- 
tion cross-section and excludes regions in the parameter space of the Minimal 
Super symmetric Standard Model (MSSM). 
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1 Introduction 

During 1998, the DELPHI experiment at LEP accumulated an integrated luminosity of 
158 pb _1 at a centre-of-mass energy, y/s, of 188.7 GeV. Results of a search for neutralino 
pair production in these data are reported here. In a separate letter [Ij], these results 
are interpreted together with those of other DELPHI searches to set mass limits on 
neutralinos, sleptons and charginos. 

In the Minimal Supersymmetric Standard Model (the MSSM) ||, there are four neu- 
tralinos Xi ) i = 1)4, numbered in order of increasing mass, and two charginos xfij = 1> 2. 
These are linear combinations of the supersymmetric (SUSY) partners of neutral and 
charged gauge and Higgs bosons. In the following, .R-parity conservation is assumed, 
implying a stable lightest supersymmetric particle (LSP) which is assumed to be the xi- 
_R-parity conservation also implies pair-production of SUSY particles, each decaying (di- 
rectly or indirectly) into a x? ; which is weakly interacting and escapes detection, giving 
a signature of missing energy and momentum. 

The neutralinos can be pair-produced at LEP2 via s-channel Z exchange or t-channel 
exchange of a scalar electron (selectron, e). The decay of heavier neutralino states to 
lighter ones typically involves emission of either a fermion-antifermion (ff) pair or a 
photon. If the scalar leptons (sleptons) are light, the two-body decay Xi $ (followed 
by I — > X°j^) mav dominate. Decays via charginos are also possible. 

Of the detectable pair production channels (i.e. excluding XiXi)-> X1X2 an d X1X3 
are important for large regions in the parameter space. For a more complete coverage, 
however, one must also consider channels like X2X3 an d X2X4) giving cascade decays with 
multiple jets or leptons in the final state. 

Moreover, a light scalar tau lepton (stau, f ) is likely to arise because of left-right 
mixing of the stau states. If the mass of the lighter stau, M fl , is close to M^o , the decay 
of the fx gives an undetectable neutralino and a low energy r which is difficult to detect. 
In this case the search for chargino pair-production has a low efficiency since xt decays 
into fxz/, but the X1X2 an d X2X2 channels are still detectable because of the r produced 
directly in the decay x\~*^\ T - It is therefore important to search also for these channels 

In the search for XkXi production with xt ~ * Xi + ff , the methods described in Refs. 
dfj were used, with minor changes. The signatures consist of pairs of jets or leptons 
with high missing energy and momentum and large acoplanarity []. In addition, several 
new searches were introduced in order to obtain a more complete coverage, in particular 
in the regions of low M^o : 

• A search for multijet events, for example from x°(i = 1, 2, j = 3, 4) with x° — > X < 2 ( ¥1 
and X2 decaying to Xi<¥l or X?7- 

• A search for multilepton events for the corresponding decays to lepton pairs. 

• A search for cascade decays with tau leptons, e.g. X2X1 production with x° — * ^ T 
and f — > X\ T 

• A search for events with low transverse energy and low multiplicity, e.g. arising from 
X2X2 production with X2 ~^ x?^ + ^~and low M^o — M^o, or from neutralino decays 
via intermediate slepton states. 

The results showed no indication of a signal and were used to derive limits within the 
MSSM scheme with universal parameters at the high mass scale typical of Grand Unified 
Theories 

1 This is defined as the the complement with respect to 180° of the angle between the jet- or lepton momenta projected 
on a plane transverse to the beam axis. 
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The DELPHI detector has been described elsewhere 0. The central tracking system 
consists of a Time Projection Chamber (TPC) and a system of silicon tracking detectors 
and drift chambers. The electromagnetic calorimeters are symmetric around the plane 
perpendicular to the beam (0=90°), with the High density Projection Chamber (HPC) 
in the barrel region (9 > 43°) and the Forward Electromagnetic Calorimeter (FEMC) 
overlapping with the Small angle Tile Calorimeter (STIC) in the forward region (1.7° < 
9 < 35°). The region of poor electromagnetic calorimetry at a polar angle close to 40° 
is instrumented by scintillators (hermeticity taggers) which serve to reject events with 
unmeasured photons. 



2 Data samples and event generators 

The total integrated luminosity collected by DELPHI during 1998 at yfs = 188.7 GeV 
was 158 pb -1 , with 153 pb _1 of adequate data quality to be used in the present searches. 

To evaluate the signal efficiencies and background contaminations, events were gener- 
ated using several different programs. All relied on JETSET 7.4 JF], tuned to LEP 1 data 
||, for quark fragmentation. 

SUSYGEN 2 . 2004 was used to generate neutralino signal events and calculate cross- 
sections and branching ratios. 

The background process e + e - -^ qq("-7) was generated with PYTHIA 5.7 JTJ. For 



/i + /i"(7) and t + t~{^), DYMU3 and KORALZ 4.2 [TT|] were used, respectively, while 



the generator of Ref. (TT2] was used for e + e — > e + e events. Four-fermion final states 



was 



were generated using EXCALIBUR [T3| and grc4f ||I4] 

Two-photon interactions giving hadronic final states were generated using TWOGAM |15 
and PHOJET while for those giving leptonic final states the generator of Ref. JT7 
used, including radiative corrections for e + e~/i + /i _ and e + e~r + r _ final states. 

The generated signal and background events were passed through the detailed simu- 
lation of the DELPHI detector || and then processed with the same reconstruction and 
analysis programs as the real data. The numbers of simulated events from different back- 
ground processes were several times the numbers in the real data, except for the number 
of simulated e + e~ pairs from two-photon interactions which was only slightly larger than 
the number expected in the data. 

In addition the simplified fast simulation program SGV, previously used in Ref. was 
adopted. SGV takes into account inefficiencies and measurement errors in the different 
tracking detectors and calorimeters, as well as multiple scattering and the showering 
of electrons and photons in the tracking volume. This made it possible to estimate 
efficiencies for points in the MSSM parameter space without full simulation, and to take 
into account all contributing production and decay channels for a given point. 



3 Event selection 

The criteria for event selection described below were based on comparisons of simulated 
signal and background event samples. The different searches used were designed to be 
mutually exclusive, in order to allow easy combination of the results. All searches used 
the information from the hermeticity taggers to reject events with photons from initial 
state radiation lost in the otherwise insensitive region at polar angles around 40° and 
140°. Events were rejected if there were active taggers in the direction of the missing 
momentum and not associated to reconstructed jets. Jets were reconstructed using the 
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LUCLUS algorithm [[| with dj oin = 10 GeV/c. Leptons were identified using the standard 
DELPHI "loose tag" criteria ||, except for electrons in the acoplanar leptons search 



(see section |3.2| ). There is a set of global event variables common to several searches. 
These were calculated based on the well-reconstructed particles in the event and include 
the total visible energy (E V - 1S ), the visible mass (M v ; s ), total momentum transverse and 
longitudinal to the beam (pt,Pl), and transverse energy (E T ). The latter is defined as 
EE", sin where Ei and 9i are the energy and polar angle of particle i. In several cases 
with two jets, their scaled acoplanarity (the acoplanarity multiplied by the sine of the 
smallest angle between a jet and the beam axis) was used. 

3.1 Acoplanar jets search 

Earlier variations of this search at lower energies have been described in Refs. @,[|. 

At least five well reconstructed charged particles were required, including at least one 
with a transverse momentum with respect to the beam above 1.5 GeV/c. The sum of 
the moduli of momenta of well reconstructed charged particles had to be greater than 
4 GeV/c, and the total transverse energy had to exceed 4 GeV. Two jets were required, 
each satisfying 10° < 6j Ct < 170° and containing at least one well reconstructed charged 
particle. Tracks which were badly reconstructed, or did not originate from the interaction 
point, were required not to carry more energy than 0A5E vis , where E vis is the visible 
energy of well-reconstructed particles. This requirement typically removes events with a 
single badly reconstructed track with a very high momentum. In addition, the calorimeter 
energy associated to such tracks had to be less than 0.2i? vis for an event to be accepted. 

Several criteria were used to reject two-photon events: the fraction of the total energy 
carried by particles emitted within 30° of the beam had to be less than 60%, the polar 
angle of the total momentum had to satisfy | cos6 p \ < 0.9, and its transverse component 
had to exceed 6 GeV/c. 

Figure |l](a) shows the distributions of invariant mass of the visible system (M v ; s ) 
divided by i/i, for real and simulated events passing the above selection. Here, and in 
the following, the simulated sample has been normalised to the integrated luminosity 
used for the data. As can be seen from the figure there is some excess of data events in 
the energy region corresponding to on-shell Z production with a lost photon from initial 
state radiation ("radiative return events"). This can be ascribed partly to a 3% deficit 
in the PYTHIA generator in this region as compared to analytical calculations ||19|| , partly 



to four-fermion processes which were not taken into account completely, and partly to 
reconstruction problems in real events with jets in the forward direction. If such excess 
events in the data pass the later steps of the selection, the background is likely to be 
underestimated and the limits derived in the absence of a signal are thus conservative. In 
the final data sample the background from Z(7) events is rather unimportant, however. 

In the continued selection events were rejected if there was a neutral particle, either 
with an energy above 60 GeV, or isolated from the nearest jet by at least 20°, and with 
an energy above 20 GeV. These criteria served to remove radiative return events. 

To reduce the WW background, events were rejected if they had a charged particle 
with momentum greater than 20 GeV/c or if the most isolated electron or muon (if 
any) had momentum greater than 10 GeV/c or was more than 20° from the nearest jet. 
Figure |I](c) shows the distributions of transverse momentum (pp) divided by i/i for real 
and simulated data, after the above selection. 

In the last step of the selection, events were accepted if they satisfied any of the 
following three sets of criteria, optimised for different neutralino mass differences (AM) . 
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The criteria involved the transverse momentum (pr), longitudinal momentum (pl), and 
invariant mass (M vis ) of the visible system, as well as the mass recoiling against it (M rec ). 
Also the acollinearity of the two jets and their scaled acoplanarity were used in this step. 
The events were accepted if: 

(i) M vis < O.ly/s/c 2 , M rec > 0.7y/s/c 2 , and p T > 7 GeV/c. In addition, the scaled 
acoplanarity was required to exceed 40°. These criteria are efficient for low AM 
(-10 GeV/c 2 ). 

(ii) 0.1^/c 2 < M vis < 0.3^/c 2 , M rec > 0.6^/c 2 , and p T > 8 GeV/c. The scaled 
acoplanarity had to exceed 25°. These criteria are efficient for intermediate AM 
(-40 GeV/c 2 ). 

(iii) 0.3v/i/c 2 < M vis < 0.5^/c 2 , M rec > 0.45^/c 2 , 12 GeV/c < p T < 35 GeV/c, and 
Pl < 35 GeV/c. The scaled acoplanarity had to exceed 25°, and the acollinearity 
had to be below 55°. These criteria are efficient for high AM (~90 GeV/c 2 ). 

Figure p] (e) shows a comparison of the scaled acoplanarity for the real and simulated 
data events passing the last step of the selection. 

3.2 Acoplanar leptons search 

The search for acoplanar leptons selects events with exactly two isolated oppositely 
charged particles (lepton candidates) with momentum above 1 GeV/c, and at most five 
charged particles in total. 

This search was slightly modified with respect to Ref. Q, as follows. The minimum 
number of TPC pad rows required for the two selected charged particles was increased 
from four to five. The lepton identification requirements were changed, accepting as elec- 
trons those particles which had an associated energy in the electromagnetic calorimeter 
exceeding half of the measured momentum, while for muons the "loose tag" criteria || 
were used. Either, both particles in the pair were required to be selected as electrons and 
not simultaneously identified as muons, or else both particles had to be muons. In ad- 
dition to the acoplanarity, the acollinearity between the two particles also had to exceed 
10°. The minimum transverse momentum required was increased from 5 to 6 GeV/c, and 
the maximum accepted energy in the STIC was reduced from 1 to 0.3 GeV. To improve 
the rejection of WW background, events with missing momentum above 45 GeV/c, and 
a scalar sum of the momenta of the two selected particles in excess of 100 GeV/c, were 
rejected. Prior to the last step of the selection 65 real data events were accepted, while 
the expected background was 62.8±4.4 events, with a contribution of 31.1±0.6 events 
from W + W~ production. As in Ref. ||, the last step involved three sets of criteria sen- 
sitive to different AM ranges. These criteria were unchanged, except for the minimal 
missing mass required in the selection optimised for large AM, which was changed from 
0.4v^/c 2 to 0.2^/i/c 2 . 

Figure |l] (b,d,f ) shows a comparison between real and simulated data for events passing 
the initial step of the above selection corresponding to rejection of Bhabha events (b), 
passing the intermediate step corresponding to rejection of two-photon events (d), and 
passing the last step (f). Real and simulated data were in good agreement throughout. 

3.3 Multijet search 

The multijet search was optimised for cascade decays of neutralinos with large mass 
splittings, giving high energy jets. Events with energetic photons, characteristic of the 
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decay x° — * Xi7> were subjected to less stringent selection criteria, giving a separate set 
of selected events with low background and comparatively high efficiency 

At least five well-reconstructed charged particles were required, and at least one of 
these had to have a transverse momentum exceeding 2.5 GeV/c. The transverse energy 
of the event had to be greater than 25 GeV, and the total energy of tracks which were 
badly reconstructed or did not originate from the interaction point was required to be 
less than 30 GeV and less than 45% of the visible energy. In addition, the calorimeter 
energy associated to such tracks had to be less than 20% of E vis . Figure |2] (a) shows the 
distributions of M vis divided by the centre-of-mass energy for real data and simulated 
background events passing the above selection. The excess of "radiative return" events 
observed in the acoplanar jets search is visible also here, and the comments of section |3J 
apply. Similarly, the deficit of events in the real data with M v ; s / y/s close to unity can be 
partly explained by a known excess of PYTHIA events with little initial state radiation. 

The total energy in the electromagnetic calorimeters had to be less than 70 GeV, 
and there had to be no single calorimeter shower above 60 GeV. The energy carried by 
particles within 30° of the beam had to be less than 60% of the visible energy. The total 
visible energy had to be less than 135 GeV, the polar angle of the total momentum had 
to satisfy | cos# p | < 0.9, and the transverse momentum had to exceed 6 GeV/c. Figure 
(c) gives a comparison of the px/v^distributions f° r re& l data and simulated background 
following the above selection. 

The scaled acoplanarity (see section |3.1| ), calculated forcing the number of jets to two, 
had to be greater than 10°. The polar angle of the most energetic jet had to be outside 
the range between 85° and 95° to avoid an insensitive detector region close to 90°, and 
its energy had to be less than 56 GeV. 

To reject WW background it was required that there be no charged particle with a 
momentum above 30 GeV/c, and that the momentum of the most isolated electron or 
muon (if any) be below 10 GeV/c, or below 4 GeV/c if the angle between the lepton and 
the nearest jet was greater than 20°. 

Events with a photon signature were then selected on the basis of reconstructed pho- 
tons in the polar angle range between 20° and 160°, isolated by more than 20° from the 
nearest charged particle track. If there was only one such photon its energy was required 
to be between 10 GeV and 40 GeV; with more than one photon, at least two had to have 
energy greater than 10 GeV. 

For the complementary sample, without a photon signature, two additional require- 
ments were imposed to reject Z7 events: the mass recoiling against the system of visible 
particles had to be greater than 100 GeV/c 2 , and all jets with energy above 20 GeV had 
to have a ratio of energy in charged particles to energy in neutral particles which was 
above 0.15. 

Lastly, events selected by the searches for acoplanar jets or leptons (sections |34] and 
|3.2| ) were rejected. Figure ^| (e) shows the acoplanarity distributions for real and simulated 
events without a photon signature passing the last step of the selection. 



3.4 Multilepton search 

The multilepton search is sensitive to cascade decays involving leptons, which can 
dominate if there are light sleptons. 

The first step in the selection, in common with the tau cascade and low searches 
(sections |3.5| and |3.6| ), was as follows. The number of charged particles was required to 
be at least two and at most eight, and events with more than four neutral particles were 
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rejected. The reconstructed invariant mass had to be below 120 GeV/c 2 , and the recoil 
mass above 20 GeV/c 2 . The calorimeter energy associated to particles which were badly 
reconstructed or did not originate at the vertex, Eb c , was required not to exceed 0.4 E vis , 
while the energy of well reconstructed charged particles had to be greater than 0.2 £vis- 
It was also required that E v[s + E\> c < 140 GeV. 

In the following step, at least two charged particles were required to be identified 
leptons. Figure |2] (b) shows a comparison between M vis / y/s distributions for real and 
simulated events passing the above selection. 

To reject Z7, two-photon, and Bhabha events, the transverse momentum of the event 
was required to exceeded 8 GeV/c, and the polar angle of the total momentum to satisfy 
I cos^pl < 0.9. The transverse energy of the event had to be greater than 25 GeV, and 
the energy in the STIC was required to be less than 10 GeV. The distributions of Pt/\^s 
for real and simulated data, following the above selection, are compared in figure (d). 

For events with exactly two isolated well-reconstructed charged particles the following 
requirements were imposed. The acoplanarity and acollinearity of these two particles had 
to exceed 15° and 6°, respectively. If the total energy in electromagnetic calorimeters 
exceeded 50 GeV the acollinearity was required to be greater than 10°. To reject W pairs 
decaying leptonically it was required that the product of charge and cosine of polar angle 
was less than —0.1 for each of the two charged particles. 

For events with two reconstructed jets, the scaled acoplanarity was required to be 
greater than 15°. 

Lastly, events selected by the searches for acoplanar jets or leptons (sections |34] and 
|3.2| ), or by the multijet search (section |3.3| ), were rejected. Figure (f) shows the distri- 
butions of acoplanarity for real and simulated data, following the above selection. 

3.5 Tau cascade search 

The tau cascade search is sensitive to X1X2 production with X2 ~^ ^ T an d f Xi T , 
where the second r produced has very low energy. The first step of the selection was 
the same as for the multilepton search (section |3.4j ), with the additional requirement of 
no more than two reconstructed jets. Two or more of the charged particles also had to 
satisfy stricter criteria on reconstruction and impact parameters. 

In the next step, the highest and second highest momenta of charged particles were 
required to be below 50 GeV/c and 25 GeV/c, respectively, and at least one charged 
particle had to have a transverse momentum above 2.5 GeV/c. Events with neutral 
showers above 300 MeV within 20° of the beam axis were rejected. The visible mass 
distributions, for real and simulated data at this stage of the selection, are compared in 
figure ^ (a). 

The criteria to reject Z7, two-photon, and Bhabha events, were the same as for the 
multilepton search (section |3.4|) , except for the minimum transverse momentum which 
was reduced to 7 GeV/c, and the removal of the transverse energy requirement. 

Figure ^| (c) shows distributions of E vis /y/s as a comparison between real and sim- 
ulated data, selected with the above criteria. There is an evident excess in the energy 
region dominated by two-photon interactions. This has been studied in a recent work- 
shop on generators at LEP2 ||20|| . The background from two-photon interactions giving 
hadronic final states is known to be underestimated, and the process 77— > £ + £~ is also 
not well described by simulation. In the case of t + t~ the treatment of tau decays in 
the generator was approximate, and polarisation effects were absent. Furthermore, some 
four-fermion processes were not completely accounted for in the simulation. If the two- 
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photon background the end of the selection is also underestimated the obtained limits 
are conservative, but in any case this background is not the dominant one. 

Events with exactly two isolated, well-reconstructed, oppositely charged particles were 
required to have acollinearity and acoplanarity above 60°. The smaller of the two mo- 
menta had to be below 70% of the greater one, and below 10 GeV/c. 

For events with two reconstructed jets the scaled acoplanarity (see section |3"H] ) was 
required to be greater than 20°, and the acoplanarity and the acollinearity greater than 
60°. 

Lastly, events selected by the searches for acoplanar leptons or jets (sections |37| and 
|3~TD or the multilepton search (section ^4|) were rejected. Figure || (e) shows the acopla- 
narity distribution for events passing the complete selection, in real data and simulated 
background. 

3.6 Low transverse energy search 

The low transverse energy (Et) search was designed to complement the multilepton 
search for cascade decays or X2X2 production with low mass splitting where X2 ~^ Xi^ + ^~ ■ 
The first step of the selection was the same as for the multilepton search. In the second 
step, it was required that there be at least three and at most five charged particles, and 
that all had momenta above 500 MeV/c. Two or more of the charged particles had to 
satisfy stricter criteria on reconstruction and impact parameters. 

In the third step, the highest and second highest momenta of charged particles were 
required to be below 50 and 25 GeV/c, respectively. At least one charged particle had to 
have a transverse momentum above 2.5 GeV/c, and at least one had to be an identified 
lepton. There had to be no neutral shower within 20° of the beam axis, and the second 
highest jet energy had to be below 30 GeV. 

Figure |3| (b) shows the distributions of / f° r events fulfilling the above criteria 
in the real and simulated data. Excess data events from two-photon interactions and the 
"radiative return" process are visible here too. Again, this could give too conservative 
limits if such excess events were to survive the complete selection. The overall effect of 
the low transverse energy search on the obtained limits is rather small, however. 

The bulk of the two-photon background was rejected by the requirements that the 
polar angle of the total momentum had to satisfy | cos# p | < 0.9, and that the transverse 
energy of the event had to be greater than 4 GeV. The distributions of M vis / yfs for the 
real and simulated data, following these requirements, are compared in figure [3] (d). 

The specific requirements for events with exactly two well reconstructed, isolated, 
charged particles were the same as in section |3.4| , with the additional requirement that 
at least one of the tracks had to have a momentum below 15 GeV/c. 

Events with transverse momentum exceeding 8 GeV / c and transverse energy greater 
than 10 GeV were rejected, unless the scaled acoplanarity, calculated forcing the number 
of jets to two, was above 20°. 

Lastly, events selected by the searches for acoplanar jets or leptons (sections |3j] and 
|3.2|) , the multilepton search (section ^4|) , or the tau cascade search (section |3.5|) were 
rejected. Figure |3] (f ) shows the distributions of scaled acoplanarity for real and simulated 
events passing the complete selection. 
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4 Selected events and expected backgrounds 

Table [I] shows the number of events selected in the different searches in real data and 
the numbers expected from the Standard Model background. Also shown are the main 
background sources contributing in each channel and the typical efficiency of each search 
for MSSM points where it is relevant. 



Search 


Data 


Total bkg. 


Main bkg. 


Typical eff. (%) 


Acoplanar jets 


19 


21.0±1.6 


W+W",ZZ 


10 - 30 


Acoplanar electrons 


16 


20.7±3.7 


W+W-,77 


10-40 


Acoplanar muons 


16 


14.6±1.3 


W+W-,77 


10-40 


Multijets, 7:s 


2 


4.3±0.5 


Z 7 


10 - 20 


Multijets, no 7:s 


39 


31.8±1.9 


Z 7 , W+W" 


10-40 


Multileptons 


23 


28.2±1.2 


W+W" 


30 - 50 


Tau cascades 


8 


9.0±1.0 


W+W",77H 


13 - 19 


Low Et 


18 


19.0±3.3 


77(— > t + t~) 


7-10 



Table 1: Results of the different searches. The typical efficiency of each search for MSSM points where it is 
relevant is shown. The efficiencies depend typically on the masses of the sparticles involved in the process. For 
any given search, events are explicitly rejected if accepted by one of the searches appearing earlier in the table. 

The main reason for the variation of the efficiencies is the variation of the masses of 
the particles involved in the process. The explicit rejection of events to avoid overlapping 
selections limits the efficiencies for those searches in which such rejection is performed 
(see section ^ and table [I]). The total number of events selected in the different searches 
was 141, with 149 ±6 background events expected. The errors given for the background 
estimates are due to the finite sizes of the simulated background samples. No error was 
assigned to account for the excesses of data events seen at early stages of the selections. 
In conclusion, the results are in good agreement with the expectation from Standard 
Model background, and no indication of a signal was found. 

5 Signal efficiencies and upper limits 

In the absence of a signal, cross- section limits were derived based on the efficiencies for 
simulated neutralino events. A total of 360 000 X1X2 events was simulated for 108 different 
combinations of masses with M^o and M^o ranging from 5 GeV / c 2 to 90 GeV / c 2 and from 
20 GeV/c 2 to 180 GeV/c 2 , respectively, and for different X2 decay modes (qqx?, 
e + e~Xi) tt ). A further 100 000 X2X34 events with cascade decays, were simulated for 
56 different points. In addition, about 5- 10 s events were simulated using SGV in order to 
obtain signal efficiencies for about 10 5 MSSM points. 

Figures [| and |5| show the expected distributions for some relevant event variables for 
X1X2 production as obtained using the full detector simulation and SGV. The efficiencies 
obtained using SGV agreed typically to ±10% relative with those obtained by full sim- 
ulation. Figure | shows a comparison between SGV efficiencies (curves) and those from 
the full simulation (points) as a function of AM in the topologies with acoplanar leptons 
and acoplanar jets. In the case of leptonic events the SGV efficiencies are generally lower, 
giving conservative limits. In the hadronic case the SGV efficiencies tend to be higher, 
and they were therefore conservatively reduced by 20% in the limit calculations. The 
effect on the M^o limit for tan j3= 1 |IJ was found to be completely negligible. 
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The limits for the X1X2 production, as obtained from the searches for acoplanar leptons 
and jets, are shown in Figs. ^ assuming different branching ratios. Similarly, Figs. |8](a,b) 
show cross-section limits for X2XA production (z= 3 or 4). For each mass combination, 
the limits were obtained by examining many possible (/x,M2) points for several tan/? 
values and high mo, where x^Xi production was kinematically allowed. The point giving 
the worst limit was taken. In the white regions marked "Not allowed", no such points 
were found. Figure 181(a) shows the limit obtained using a Bayesian combination ETI] of 



the results from the multijet and acoplanar jet searches in the case where x°i~ an d 
X^^XiQQ- Figure §(b) gives the corresponding limits when xH~ *x\li as obtained from 
the search for multijet events with a photon signature. 

In addition to such limits on the production cross-sections, the approach using a fast 
simulation makes it possible to scan regions of the MSSM parameter space and calculate 
the efficiencies directly at each point, simulating all neutralino production channels and 
decay chains. Since they were defined to be mutually exclusive, the different selections 
can be combined using the Bayesian multi-channel approach [21] to obtain the exclusion 
confidence level for each set of MSSM parameters^. Figs. |9] and [K] show the regions 
excluded by the different contributing searches in the (/x,M 2 ) plane for tan (3 = 1 and m = 
ITeV/c 2 and 80 GeV/c 2 , respectively. Also shown are the combined exclusion regions 
for the two values of tuq. In the region indicated as "Not allowed" the lightest chargino 
is lighter than x\- Although the process for which it was designed is not important here, 
the r cascade search is efficient for cascade decays involving leptons in the region close to 
fi = for itiq = 80GeV/c 2 , and when the chargino- neutralino mass difference is small. 
(In the latter case the decay X23 — * Xi^ u 1S followed by an almost invisible chargino 
decay.) 

The thin dotted curve in the figures indicates the chargino isomass contour corre- 
sponding to the kinematic limit for X1X1 production. For high m this is very close to 
the exclusion limit from chargino searches. For low mo the region excluded from chargino 
production is smaller |[|, but the neutralino excluded region is increased, as can be seen 
from Figure [H| Therefore the overall limit on M^o for tan/5 = 1 is determined by the 
intersection of the chargino isomass contour with the region excluded by neutralinos for 
high mo 0. The corresponding Xi isomass contour is shown as the thin dashed curve. 

At low itlq and low M 2 , the region excluded by neutralinos shrinks with increasing 
tan (3 due to enhancement of the invisible x\ & ~^ u Xi decay channel. There is no 

substantial change of the high m exclusion region with the increase of tan j3. 



6 Summary 

Searches for neutralinos at y/s = 188.7 GeV, using several mutually exclusive sets of 
criteria, gave no indications of a signal. As a consequence, upper limits on cross-sections 
for different topologies were derived, ranging from about 0.1 pb to several picobarn. The 
efficiencies computed with a full simulation of the DELPHI detector were extended to 
the whole range of the SUSY parameters explored by using a fast detector simulation, 
which included all neutralino production and decay channels. Exclusion regions in the 
MSSM parameter space were then derived. The methods used were designed for deriving 
general MSSM mass limits in the Minimal Supersymmetric Standard Model, as done in 
a separate letter M . 

2 The same procedure is applied in Rcf. Ill, including also the production of other supersymmetric particles. 
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Figure 1: The comparison between the real and simulated data for the acoplanar jet 
selection (a,c,e) and acoplanar lepton selection (b,d,f) is shown. Plots (a,b) show the 
visible mass divided by the centre-of-mass energy at an initial stage of the selections. 
Plots (c,e) shows the missing transverse momentum divided by centre-of-mass energy at 
an intermediate stage of the selections. Plots (e,f) show acoplanarity distributions after 



the last step of the selections. The selections are described in sections |3J] and |Q| Plots 
(e,f) also show the expected signal of X1X2 production for two different neutralino mass 
combinations assuming a cross-section of 1 pb and the decay X2 ~^ %*Xv 
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Figure 2: The comparison between the real and simulated data for the multijet selection 
(a,c,e) and multilepton selection (b,d,f ) is shown at three different stages of the selection. 
Plots (a,b) show the M vis divided by the centre-of-mass energy at an initial stage of the 
selections. Plots (c,d) show the missing transverse momentum at an intermediate stage 
of the selections. Plots (e,f) show the acoplanarity after the last step of the selections. 
The selections are described in sections [373] and |37J]. The distributions expected for X2X3 
production with X3^xffi~ ¥ Xi^' 1 normalised to a cross-section of 2 pb, are also shown 
for decays into quark and lepton pairs in e) and f), respectively (dashed histograms). 
Equal mass differences — .!/ •.• M^o — M^o 25 GeV/c 2 were assumed. 
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Figure 3: The comparison between the real and simulated data for the tau cascade 
selection (a,c,e) and low E T selection (b,d,f) is shown. Plots (a,b) show the M vis divided 
by the centre-of-mass energy and transverse momentum divided by the centre-of-mass 
energy at an initial stage of the selections. Plots (c,d) show the visible energy divided 
by the centre-of-mass energy at an intermediate stage of the selections. Plots (e,f) show 
the acoplanarity after the last step of the selections. The selections are described in 
sections |3.5| and |3.6| . The dashed line in e) shows the tau cascade signal expected from 
XiX°2 production with xl -» rf -> ttx\ , A% = 34.8 GeV/c 2 , M f = 36.8 GeV/c 2 , 
and M^o = 60 GeV/c 2 . In f) the dashed line corresponds to X2X2 production with 
= e,^,r), M^o = 35 GeV/c 2 , and M^o = 40 GeV/c 2 . The signals are 
normalised to 2 pb. 
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Figure 4: The expected distributions of relevant event variables characterizing the X1X2 
production with M^o = 40 GeV/c 2 and M^o = 80 GeV/c 2 , as obtained using the full 
detector simulation (DELSIM) and SGV for the acoplanar jet (jj) topology (upper four 
plots) and acoplanar lepton (11) topology (lower four plots). The decays X2 ~> XiQQ or 
X®£ + £~ were assumed as appropriate (£ + £~ denotes e + e~ and in equal proportions). 
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Figure 5: The expected distributions of relevant event variables characterizing neutralino 
production for the multijet topologies (upper four plots) and the tau cascade topology 
(lower four plots). In the multijet case chosen, X4X2 production dominates with 50% of 
the x° decaying to Xil- The neutralino masses are M^o = 31 GeV/c 2 , M^o = 60 GeV/c 2 , 
and M^o = 100 GeV/c 2 . In the tau cascade case, X1X2 production with X2 ~^ T ^ —* TT Xi 
was assumed with M^o = 34.8 GeV/c 2 , M f = 36.8 GeV/c 2 , and M^o = 60 GeV/c 2 . 
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Figure 6: Efficiencies for XiXi production as obtained with the full simulation (DELSIM, 
points with error bars) and the fast simulation (SGV, points connected by straight lines) 
for different AM — M^o—M^o, assuming the decays x°— >x°ff (f=/i,e,q). The mass of xl 
was fixed to 35 GeV/c 2 . 
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Figure 7: Contour plots of upper limits on the cross-sections at the 95% confidence 
level for X1X2 production at y/s = 189 GeV. In each plot, the different shadings 
correspond to regions where the cross-section limit in picobarns is below the indicated 
number. For figures a), b), c), X2 decays into Xi an d a) e + e~, b) and c) qq, while 

in d) the branching ratios of the Z was assumed, including invisible states. The dotted 
lines indicate the kinematic limit and the defining relation M^o > M^o. 
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Figure 8: Upper limits on the cross-sections at the 95% confidence level for x^Xi pro- 
duction with x® — tX^QQ (^=3,4) at y/s = 189 GeV. The different shades correspond to 
regions where the cross-section limit in picobarns is below the indicated number. In the 
darkest shaded regions there are points which are not excluded for any cross-section. x\ 
was assumed to decay into XiQQ i n a ); an d into Xil i n b). The limits in a) are based on 
the acoplanar jets and multijets selections, while those in b) derive from the search for 
multiiets with photons. 
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Figure 9: Regions in the (//, M 2 ) plane excluded at 95% confidence level for tan (3—1, 
assuming mo = 1 TeV/c 2 . The exclusion by individual searches for jet pairs (top left), 
multijets without 7:s (top right), multijets with 7:s (bottom left), and r cascades (bottom 
right) are compared with the combined exclusion based on all searches (thick dashed 
curve and light shading). Also shown are the kinematic limit for chargino production 
(thin dotted curve) and the isomass contour for the minimum allowed neutralino mass 
]Ij (thin dashed curve). In the region marked "Not allowed" the chargino is the LSP. 
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Figure 10: As figure but for tuq = 80 GeV/c 2 and six different contributing searches. 
From left to right and top to bottom: jet pairs and multijets without 7:s (hatched), 
multijets with j:s, lepton pairs and multileptons (hatched), and r cascades. 



